The Bose-Fermi Kondo model with a singular dissipative spectrum: Exact solutions 

and their implications 
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Quantum dissipation induces a critical destruction of a Kondo screened state, which is of interest 
in the contexts of quantum critical heavy fermion metals and magnetic nanostructures. The sub- 
ohmic Bose-Fermi Kondo model provides a setting to study this effect. We find that this many-body 
problem is exactly solvable when the spectrum of the dissipative bosonic bath, J{ijj), is singular, 
corresponding to J(r) = const.. We determine the local spin correlation functions, showing that the 
singular longitudinal fluctuations of the bosonic bath dominate over the transverse ones. Our results 
provide evidence that the local quantum critical solution, derived within the extended dynamical 
mean field approach to the Kondo lattice model, has a zero residual entropy. 
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The Bose-Fermi Kondo model (BFKM)[1 i has be- 
come of considerable interest in a number of different 
contexts, including the extended dynamical mean field 
approach to the heavy fermion quantum criticalit\_ 
the dissipative effects in mesoscopic structures 0,15] 
model describes a quantum impurity spin coupled simul- 
taneously to a conduction-electron bathQ and to a dis- 
sipative bosonic continuum[3, S|- The coupling of the lo- 
cal moment to conduction electrons leads to the Kondo 
screening effect. The coupling to the bosonic bath, on the 
other hand, causes the critical suppression of the Kondo 
effect. The latter is the key ingredient of the local quan- 
tum critical solution to the Kondo lattice model0], and 
has received some fairly direct support from the exper- 
iments in heavy fermions[^, [l^. In even broader con- 
texts, models such as BFKM serve as prototype systems 
to study quantum dissipation and quantum coherence 

The BFKM has a rich phase diagram, which is sensitive 
to the form of the spectral density of the bosonic bath: 



J{u) 



sgn(w) 



(1) 



Most of the previous studies have focused on the sub- 
ohmic case with < e < 1. The perturbative renor- 
malization group (RG) analysis [l|, 0, [ll|, 12 1 shows 
that spin-anisotropy is a relevant perturbation, and in- 
deed the Kondo-destroyed phase of the Ising-anisotropic 
BFKMP, [m, [H, El, m, [3 differs qualitatively from 
that of the spin-isotropic BFKM [13, [ill. The cases of 
singular bosonic spectrum [J{lo) being divergent in the 
zero- frequency limit], e > 1, have received much less 
attention. (A singular bosonic spectrum might be re- 
alizable in some low dimensional ferromagnetic single- 
electron structures [^.) One may argue [l^l on general 
grounds that the large bosonic spectral weight at low 
energies implies that an infinitesimal coupling between 
the local moment and the bosonic bath will suppress the 
Kondo effect. For the Ising-anisotropic BFKM, the nu- 
merical renormalization group (NRG) work of Glossop 



and Ingersent [l6| has indeed demonstrated this. 

Several factors have motivated us to search for an exact 
solution to BFKM. It will put the above perturbative, nu- 
merical, or saddle-point results on a more firm ground, 
thereby elucidating the physics of not only the Kondo 
destruction but also the quantum dissipation/coherence 
in general. In addition, it may help resolve some ap- 
parent inconsistency concerning the entropy of the local 
quantum critical solution to the Kondo lattice system. 
Some general considerations suggest a zero residual en- 
tropy (see below), yet Kircan and Vojta [11] have inter- 
preted their results for a dynamical large-N limit of the 
spin-isotropic BFKM as signaling a contrary conclusion. 

In this letter, we show that the BFKM with a sin- 
gular bosonic bath, corresponding to e = 2, is exactly 
solvable. The exact solution establishes that an infinites- 
imal coupling between the local moment and the bosonic 
bath indeed causes a Kondo destruction. Surprisingly, we 
find that the spin-isotropic problem has the same univer- 
sal properties as the Ising-anisotropic one, which implies 
that the longitudinal dissipative coupling plays a dom- 
inant role. This effect turns out to be missed by the 
large-N limit of Ref. [lO] in an interesting way. 

The Hamiltonian for the BFKM is 



H — Hp + Hb + Hj + Hq + Hh- 



(2) 



■ £kc\^Cka and Hb 



Here, Hp — J2k,a ^k^ka^^a "i^^ ^^b — L^p^a "^P^pa-^pa 

respectively describe a conduction-electron band and a 
bosonic bath, and and Wp their corresponding disper- 
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{Jj2)[s-{0)S-^ 



-iO)S-] + J||s^(0)5^ 



specifies the Kondo coupling between a spin- 1/2 local 
moment, S, and the conduction electron spin Sc = 



E I'^fecr^cr.o-'Cfccr' (wlth T bclug the Pauli matrices) and 
Hg = Qa^'^S'^ the impurity coupfing to the bosonic 

field, with = Yjpi^pa + 'i'Lpa) at the impurity site 
0. The dissipative bosonic spectrum, J{^JJ) = ^p[5{ijJ — 
Wp) — 5{u> + Wp)], is taken to have the form of Eq. ([T|). 
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In order to probe the local spin responses, introduces a 
local magnetic field, h, via Hh = hS^. 

The partition function can be written in the path- 
integral form. The quantum impurity spin S is repre- 
sented by an SU[2) coherent state S^(t) with the con- 
straint r2^(r) ~ 1. A Berry phase term iSb[J^] character- 
izes the quantum dynamics of the impurity spin. Mean- 
while, the bosonic bath can be traced out, leading to a 
long-ranged interaction (along the r-dimension) for the 
impurity spin [3, Q: 



(3) 
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X exp{^ / / 
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dTldT2 
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with Zb being the partition function of the bosonic bath. 
(Note that the purely classical model defined along a 
chain, < t < /?, is considered to be ill-defined for 
1 < e < 2, since the long-ranged interaction makes the 
ground state energy per unit length diverge in the /3 — > oo 
limit, but the quantum problem we consider is still well- 
defined. [2lj) For e = 2, the long-range interaction in 
the last term becomes Yla=i[Io ^'''9a.S^''iT)/2]'^, "^hich 
can be further decomposed by introducing a Hubbard- 
Stratanovich vector A = (Ai, A2, A3): 



exp{^[ / dTSgan^{T)/2\'} = 

a=l -^0 

r. 3 .0 

7r-3/2 / rfXcxp{-A2-^g„A„ / dr51)"(T)}. (4) 

The total partition function of the BFKM, when rewrit- 
ten in the operator formalism, becomes 



dAe-^ Tre-'^^^^), 



U(\) = Hp + Hj + Hh 



(5) 



where ~ J2a=i 9a^aS°'. Hence, we are led to solve the 
BFKM, Eq. in terms of a pure Fcrmi-Kondo model 
in the presence of a Gaussian-distributed magnetic field 
A, along with the external field h. 

The Bose-Kondo model. We first consider the case 
with Qa = g (for a = 1, 2, 3) and in the absence of con- 
duction electrons, i.e., the SU{2) Bose-Kondo model. 
Here, H{X) = hS^ + gX ■ S. Using the parameter- 
ization A = (A sin(p COS0, A sin (y9 sin0, A cos we can 
specify the eigenvalues of H{X) as ztgX/2, with A ~ 
A-^(l — y)2 + 4ycos2 ^ and U — jx- Inserting these 
into Eq. ([5]) leads to the impurity partition function, 

Zloc 
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2 cosh 



hp g^P sinh ■ 



(6) 



It follows that the static local spin susceptibility is 



Xlo 



1 dHnZia 
(3 dh^ 



h=0 



12 1 + 52/52/8 
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(7) 



The asymptotic behavior in the low-temperature limit 
has a Curie form, xioc{T — > 0) = /3/12, with a reduced 
Curie constant (1/12 instead of the free-spin value, 1/4). 

The dynamical local spin-spin correlation function, 
XIoc{t) = {S^{t)S^{0))ioc, is obtained from a Gaus- 
sian averaging: xioc{t) = (n^^^^/Zioc) J dXe^^'A{X), 

where A{X) = Tre-'3^[^lS'^(r)S'^(0). The trace 
amounts to E„.,„ e-^"''e(^''-^™)"|(n|^^|m)|2, 



where 
states 

|-> = 
A{X) - 
turn, 



n 
of 



and 
H[X]-- 
(- sin 



i[cos2<^cosh^ 



over all the eigen- 
= (cos |^e~*2 ^ sin ^e*2 ), 

It follows that 

gA(/3-2T) 



, cos f e*2 ). 
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In the asymptotic low-temperature and long-time limit 
(r 13/2, (3 — > 00), it also has a Curie form. The full 
result is also illustrated in Fig. [TJ 
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FIG. 1: The dynamical loci spin susceptibility of the SU{2) 
Bose-only Kondo model at e = 2 (with g — 1 and (3 = 10). 

We next consider the Ising case, i.e., ffi = 52 = 0, gs = 
g. The impurity partition function is now 

= 2cosh^e^. (9) 

Both the static and dynamical local spin susceptibility 
have simple Curie forms: xioc = /3/4, and xiociT) = 1/4. 

To summarize, the local spin response of both the 
SU{2) and the Ising cases has a Curie form in the low 
temperature hmit. For the SU{2) problem, the fact that 
it has the same universal behavior as the Ising case im- 
plies that it too is controlled by the g* = 00 fixed point; 
in other words, the longitudinal fluctuations dominate 
over the transverse ones. This is in contrast to the e < 1 
case, where the SU{2) and Ising problems are controlled 
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by a finite-q* fixed point and a g* = oo one. respectively. 

S El, 121 Our exact result for the e = 2 case reveals 
the link between a singular bosonic spectrum and domi- 
nant longitudinal fluctuations. We therefore expect that, 
for the 1 < e < 2 cases too, the SU{2) problem will be 
controlled by the fixed point at g* = oo. 

The Bose-Fermi Kondo model. We now incorporate 
the Kondo coupling as well. The spin-isotropic case re- 
quires the usage of Bethe-ansatz method, which will be 
discussed elsewhere [22|. Here we consider the case in 
which the pure fermionic Kondo part is placed at its 
Toulouse point and, moreover, the impurity-boson cou- 
pling is purely Ising. Since the spin-anisotropy in the 
impurity-bosonic coupling is irrelevant and our results 
will demonstrate that an infinitesimal coupling to the 
bosonic bath leads to the Kondo destruction, the spin- 
isotropic BFKM should behave similarly. 

When the Kondo couplings take the Toulouse val- 
ues, the model Eq. (O can be mapped^ to the follow- 
ing non-interacting spinless resonant level model (RLM): 



1 r"' 
Trr r^+pA^ 



4 arctan^ ^ 



Ht{\) = Efc 



+ V{cU + h.c.) +ed{dU-\) 



Here, the impurity spin is represented by S*^ d^ d—l/2, 
d^ , and S~ — > d, with d describing a spinless 
fcrmion. The hybridization V is proportional to J_l, 
while Ed ~ gX + h. After integrating out the fermionic 
bath, we obtain Zt = ir"^/^ J dXe'^^Tie-'^"^^^'' = 
Zc^ioc- Here, Zj, is the partition function of the elec- 
tron bath, and Zioc = 



Z,o,[A] with 



— n{Lu) arctan } 



•cosh(— ), 



(9) 



where n(w) = [1 -I- exp(/3<x')]~^ is the Fermi function, F = 
i^PoV"^ the bare resonance width, and po the conduction- 
electron density of states at the Fermi energy. Without 
loss of generality, we assume a flat band for the electrons, 
taking the usual limit of a large bandwidth D F). 

Eq. ^ has been derived with care, containing not only 
the usual phase shift contribution [y], but also a resid- 
ual atomic term. Indeed, it recovers the results for both 
the Ising bosonic Kondo model (F = 0) and the conven- 
tional RLM {g = 0). When both couplings F and g are 
non-zero, the Gaussian averaging over A complicates the 
problem. We focus on the zero-field static local suscepti- 
bility. It is straightforward to show that 



Xlo 



/_°°^dAe-^^Xioc[A]Z;„JA], 
/r^dAe->^Z,„,[A] 



(10) 



where xioc[X] - { J^A/joc [A] + /JA/^JA]} and ^^.[A] = 
^OhiZip|A]^^ At low temperatures, we use the asymptotic 
approximation n{uj) ~ @{—uj) (Sommerfeld expansion, 
valid up to corrections of order of T^/F) and integrate 



over uj, obtaining Xioc[A] 

and Zioc[\] = exp{^ arctan^ ~ + {^?]}Z', 

where Z' is independent of A in the limit D /T ^ oo. 

In the weak boson-impurity coupling regime, g/T <^ 1, 
the A- integrals in Eq. fTUIl can be asymptotically ap- 
proximated by introducing a cut-ofF A ?» T/g ^ 1. 
Within the cut-ofF (|A| <C A, called region I below), 
the dummy variable ^ is always small such that ex- 
pansions like arctan ^ ^ are appropriate. Outside 
the cut-ofF (|A| ^ A, region II), ^ is large enough and 
one has arctan^ « sgn(A)^. Xioc[A] then becomes 

and I For |A| < A and |A| > A, re- 
spectively. .Zioc [A] contains a dimensionless combination, 

a 2 2 

a = A characteristic temperature scale, T* = 

(corresponding to a = 1), arises, separating two limiting 
temperature regimes with distinct asymptotic behavior 
for Xioc- In the low temperature regime, i.e., T <^T* ( 
or a 1), the A-integration is dominated by the region- 
II contribution, and x/oc ~ -f. In the higher temperature 
regime, T* < T < F (or a ^ 0), on the other hand, 
the rcgion-I contribution dominates, and xioc ~ A 
smooth, but broad, crossover takes place around T ^ T* . 
The phase diagram is illustrated in Fig. [S] 

In short. For a fixed non-zero Kondo coupling (and, 
hence, non-zero F), the local spin susceptibility in the 
low-temperature limit turns to a Curie form for any non- 
zero g. This implies that an infinitesimal g causes a de- 
struction of the Kondo effect; the BFKM is controlled 
by the fixed point associated with the Bose-only Kondo 
model solved in the previous section. 
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FIG. 2: The phase diagram of the present BFKM. xioc has 
the Curie or Pauli forms when T < T* or T* < T < F; for 

«r, r* «ffV(27rr). 

Implications for the EDMFT solution of the Kondo 
lattice model. The extended dynamical mean field ap- 
proach to the Kondo lattice model is studied through 
a BFKM with self-consistent spectra for the bosonic 
and fermionic baths. In the local quantum critical 
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solution,!^, [l^, 23, 24 1 the self-consistent bosonic bath 
has a spectrum with e = 1^. 

For e < 1, we have the standard sub-ohmic situation in 
which an unstable fixed point (critical, referred to as "C" ) 
separates the Kondo and LM fixed points. [13] NRG 
studies of the Ising- anisotropic BFKM[i3| found that, for 
1 < e < 2, an infinitesimal g leads to Kondo destruction. 
(Similar results were inferred from an 0(N) representa- 
tion of the Ising-anisotropic BFKM[i^.) In other words, 
at e exactly equals to 1, "C" and the Kondo fixed points 
are merged and the impurity entropy will vanish. Wc 
expect the impurity entropy to be a smooth function of 
e (which is consistent with the NRG result [i^), so the 
impurity entropy should vanish at e = 1~ as well. It fol- 
lows that the corresponding solution to the Kondo lattice 
has a vanishing residual entropy, which has been experi- 
mentally seen in YbRh2Si2 , a prototype quantum critical 
heavy fermion metal. [26| . 

The situation is very different in the large-N limit dis- 
cussed in Ref. [2^ , where "C" and the Kondo fixed points 
are still separated even for 1 < e < 2.[23] Moreover, 
the impurity entropy, at both the "C" and the LM fixed 
points, are finite for the entire range of < e < 2, in- 
cluding at e = 1. 

Our exact results imply that the large-N limit becomes 
problematic when the bosonic bath is singular. The lo- 
cal susceptibilities of both the LM and "C" fixed points 
of the large-N limit were shown to have xioc{t) ~ 1/t'', 
with r/ = e for the entire range of < e < 2 the 
structure of fixed points remains the same over this range 
< e < 2 This would give rise to ry = 2 for e = 2. 
But our exact solution for the SU{2) case [c/., Eq. ([5]) 
and Fig. [T] shows that 77 = for e = 2. The problem lies 
in the large-N- limit's under-estimation of the longitudi- 
nal fluctuations, which is of order 1/7V. As we discussed 
earlier, when the bosonic bath has a singular spectrum 
(e > 1) and for any finite-N, the strong longitudinal fluc- 
tuations dominates over the transverse fluctuations, and 
the spin-isotropic problem fall in the same universality 
class as the Ising case. This implies that the structure of 
the fixed points of the spin-isotropic BFKM changes as e 
passes through 1. 

Keeping track of the longitudinal fluctuations is also 
important to enforce the Griffiths' bound 28|, which re- 
quires XIoc{t) not to decay faster than the interaction, 
l/r^"*^. This bound is violated (for e > 1) by the N = 00 
result. [20j xiociT) ~ A contribution of order l/r^"*^ 

to xiociT) occurs at order 1/iV, which is missed in the 
N = 00 limit. On the other hand, our exact result for 
the SU{2) problem at e = 2, XIoc{t) ~ satisfies the 

Griffiths' bound. 

In conclusion, we have presented the first exact solu- 
tion to the Bose- Fermi Kondo model. Our results display 
the physics of dissipation-induced Kondo destruction, an 
effect that is of extensive current interest in the under- 
standing of quantum criticality in heavy fermion metals. 



Our results also establish the important role that the 
longitudinal fluctuations of the dissipative bath plays in 
the Kondo destruction, which has been missed in an in- 
teresting way in a large-N formulation of this problem. 
The new understandings achieved here have allowed us 
to elucidate the issue of zero residual entropy of the local 
quantum critical solution to the Kondo lattice models. 

We thank R. Bulla, M. T. Glossop and, particularly, 
M. Vojta, for useful discussions, and the NSF Grant No. 
DMR-0706625 and the Robert A. Welch Foundation for 
support. One of us (J.D.) has also been supported by the 
NSF of China and the PCSIRT(IRT-0754) of the Chinese 
Ministry of Education. 



[1] 
[2: 

[4: 
[6: 

[7 



[9 

[10: 

[11 
[12: 

[13: 

[14 

[is: 
[16: 

[17 

[is: 

[19 
[20 

[21 



[22] 
[23] 

[24] 
[25] 
[26] 
[27] 



J.L. Smitli and Q. Si, Europhys. Lett. 45, 228(1999); Q. 

Si and J.L. Smitfi, Pliys. Rev. Lett. 77, 3391(1996). 

A.M. Sengupta, Phy. Rev. B 61, 4041(2000). 

Q. Si, S. Rabello, K. Ingersent, and J.L. Smith, Nature 

(London) 413, 804(2001); Phys. Rev. B68, 115103(2003). 

K.Le Hur, Phy. Rev. Lett. 92, 196084(2004). 

S. Kirchner et ai, PNAS USA 102, 18824(2005). 

A.C. Hewson, The Kondo problem to Heavy Fermions, 

(Cambridge University Press), United Kingdom, 1993. 

A.J. Leggett et al, Rev. Mod. Phys.59, 1(1987). 

S. Sachdev and J. Ye, Phys. Rev. Lett. 70, 3339(1993); 

A. J. Bray and M. A. Moore, J. Phys. C13, L655 (1980). 

H. von Lohneysen, A. Rosch, M. Vojta, and P. Woffle, 

Rev. Mod. Phys. 79, 1015(2007). 

P. Gegenwart, Q. Si, and F. Steglich. larXiv: Q712.2045] 
L. Zhu and Q. Si, Phy. Rev. B 66, 024426(2002). 
G. Zarand and E. Dernier, Phy. Rev. B 66, 024427(2002). 
D.R. Grempef and Q. Si, Phys. Rev. Lett. 91, 
026401(2003); J.-X. Zhu et al, ihid. 91, 156404(2003). 
R. Bulfa, N.H. Tong, and M. Vojta, Phys. Rev. Lett. 91, 
170601(2003); M. Vojta et al, ibid. 94, 070604 (2005). 
M.T. Glossop and K. Ingersent, Phys. Rev. Lett. 95, 
067202(2005) 

M.T. Glossop and K. Ingersent, Pliys. Rev. B 75, 
104410(2007). 

M. Vojta, C. Buragohain and S. Sachdev, Phys. Rev. B 
61, 15152 (2000). 

L. Zhu et al. Phys. Rev. Lett. 93, 267201(2004). 
Q. Si, unpubfished notes (2004). 

M. Kircan and M. Vojta, Phys. Rev. B 69,174421(2004). 
Note however the overalf pre-factor ^ in the parti- 
tion functions, Eqs. (|6I9|I . Its origin fies in the renormaf- 
ization of by g. Some regufarization procedure may be 
introduced to cancel this pre-factor. This, however, does 
not affect the focal spin responses. 
J. Dai, C.J. Bofech, and Q. Si, unpubfished (2007). 
M.T. Gfossop and K. Ingersent, Phys. Rev. Lett. 99, 
227203(2007). 

J.X. Zhu et al, Phys. Rev. Lett. 99, 227204(2007). 

R. Buffa, private communication (2007). 

J. Custers et al, Nature 424, 524 (2003). 

The refevant farge-N cafcufations for the Bose-Kondo 

probfem (J = 0) were reported in Ref. [lO]. Those for 

the BFKM (finite J, incfuding at the unstabfe "C" fixed 



5 



point) were done (using the formalsim of Ref. (3) in L. [28] R. B. Griffiths, J. Math. Phys. 8, 478 (1970). 
Zhu and Q. Si, unpublished notes (2004). 



